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Abstract 

We present an algorithm to find an integral vector in the polyhedral cone T = {X\AX < 0}, 
without assuming the explicit knowledge of A. About the polyhedral cone, V, it is only 
given that, (i) the elements of A are in {— d, — d+ 1, • ■ • , 0, • ■ • , d — 1, d} r d e N, and, (ii) 
Y = [y(l),y(2), ■ ■ ■ ,y(n)] is a non-zero integral solution to T. The proposed algorithm finds a 
non-zero integral vector in T such that its maximum element is less than (2d) 2 " ~ 1 /2™~ 1 . 

Index Terms — Integer programming, Algorithm. 



I. Introduction 

Finding an integral vectoi0 in the polyhedral cone T = {X\ AX < 0}, for a given/known 
matrix A, A G Z nxn , is a problem which has been considered in great detail (U, @, 

HI. 

In this paper, we consider the above-mentioned problem from another angle with a 
distinctly different assumption. Here, we assume there is no explicit knowledge about 
A (i.e., A is unknown), but a non-zero integer solution of T is given. Under these 
assumptions, we can show that not only does there exist another integer solution of 

T but also the maximum element of the obtained solution is less than - ' 2n _ 1 — , when 
elements of A are in {— d, —d + 1, ■ • • , 0, • ■ ■ , d — 1, <i}H, where d G N. 

The rest of the paper is organized as follows. In Section [TT] we present the main 
theorem of the paper. In Section [III] we consider an example. In Section [TV] we present 
some applications of the algorithm. 

II. Main Theorem 

In this section, we present the main theorem of this paper. The proof of this theorem 
is constructive and along the proof we present the algorithm that archives the desired 
properties. 

1 An integral vector is a vector with non-negative integer elements. 

2 In other words, suppose that Y — [y(l),y(2), . . . ,y(n)] is a given non-zero integer solution of V such that 
maxi<i<„ y{i) > ^ 2d ^„_ 1 — . Then, without knowing A, we can find another non-zero integer solution of Y such 
that its maximum element is less than ^ 2d \,,_ , . 



November 15, 2011 



DRAFT 



Theorem 1: Consider Y = {X|AX < 0}, with this knowledge that the elements of 

A are in {— d, — d + 1, • • • , 0, • • • , d — 1, d}. Assume that Y = [y(l), y(2), . . . , y{n)\ is a 
given (arbitrary) non-zero integer solution of Y. Then there exists an integer solution 
X = [x(l),x(2), . . . ,x(n)j which satisfies AX < and 



maxa;(i)< — 7 ^ zl — . (1) 



(2d) 2 

Proof: Without loss of generality we assume that, 

y(l)<y(2)<...<y(n). 

We consider the polyhedral cone Ay, Ay c Y, which is defined by all inequalities of the 
form 



^ Cjxji) < 0, 



1=1 



where G {—d, ■ ■ ■ ,d} for all % and YTi=\ C %V^S) ^ 0- 

We will describe a procedure to construct a solution X G Ay that satisfies 



(2d) 



2" 



max < 



n-l 



Definition 2: The sequence Tj is defined by the recurrence relation Ti — d, T 3 ■. = 
2T 2 
Thus 

T, = \{2df-\ 

For j G {1, . . . , n — 1}, let A y denote the polyhedral cone in n+ 1 — j dimensions which 
is defined by all constraints of the form 

n 

^2 < 0, 

i=j 

where for all i, q is an integer with |q| < T jr and Yl7=j c iV(^) — 0- 

Our procedure begins by setting x(n) = 1. For integers j decreasing from n—\ down to 

1, we describe a way to select and update the partial solutions, that is, in each iteration, 

say n — j, we derive a feasible solution with lower dimension called partial solution in 

that iteration and denote by X^ n ~^= {x(i) : % > j} G A y . 

For j — n — 1 and x(n) = 1, any real partial solutions X^ 1 ' 1 G A y _1 is calculated by 

constraints of the form 

tf-ix(n-l) < ^ (2) 

or 

- 1) > V n , (3) 
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where 7™_ 1 , 7n> lli-n an< ^ In are integers satisfying 

< 7n-l < T n _! (4) 

1 < 7^ < (5) 
1 < 7i_i < T n _! (6) 

< 7n < T n _! (7) 

7n-i^-l) < 7^(n) (8) 
j n -iV(n-\) > ~i l n y(n) (9) 

We initially choose x(n — l) to be a positive number that satisfies all constraints in A y 1 
assuming that x{n) = 1. More specifically we choose positive integers 7*_ x < T n _i and 
7^ < T„_i such that 

7* 

J 1 = x(n — 1) < x(n) = 1, 

7n— 1 

and the constraints of A y _1 are satisfied. We next multiply x(n — 1) and x(n) by 7*_ x to 
obtain an integral partial solution satisfying 



1 <x(n-l) < x(n) < 7*_! < T 



n— 1) 



and this will be our initial partial solution when we begin to consider x(n — 2). 

For j decreasing from n — 2 down to 1, suppose that we have an integral partial 
solution X("-^ +1 ))* = : i > j + 1} e A y +1 with x*(n) < FEj-i T *- We wil1 use this 

partial solution to construct an integral partial solution X^ n ~^ e A y with < 11?=/ 
We begin by setting = j + 1 < i < n. Assuming that this is a legitimate 

assignment, in order for real partial solution X^ n ~^ to be an element of A y , we must 
have that 

n 

^ Cjxji) < 0, 

i=j 

where for all i, q is an integer with |q| < Tj, and Y17=j C «2/W ^ 0- There are three cases 
to consider for cf 

1) If Cj = 0, then since X^ +1 »* e A y +1 it follows that YILj+i c i x( S) < holds when 
Ci is an integer with |q| < Tj for all % > j + 1, and when X)" =J - +1 c;y(i) < 0. 

2) If cj = dj > 0, then we obtain an upper bound on x(j): 

1 n 

Observe that 



°J i=j+i 



y(j) < ~vu E ^(0- (id 



c 
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3) If Cj = c l j < 0, then we obtain a lower bound on x(j): 

1 n 

X U) > -T E (12) 

Observe that 

2/0') > -4 E (1 3 ) 

In order for this approach to lead to a valid partial solution, we need to guarantee that 
all upper bounds on x(j) exceed all lower bounds on x(j). By (flCfy - pS]) we want to 
establish that 

i=j+l 3 i=j+l 

when 



Constraint (|14|) is equivalent to the condition 

n 

E(4i-c?4K«<0. (16) 
The property (|T5|) can be rewritten 

E ( c ' c ? - ^ °- (1 7 ) 



Notice that since \c\\ < Tj and |c"| < T 3 - for all z > j, it follows that |qcj — c"cj-| < 2Yj = 
T i+ i for all % > j + 1. Since holds and X( n ^' +1 »* e A^ +1 by assumption, it follows 
that (dU) holds. 

We choose to be the maximum value satisfying all constraints (|T0|) and ((12)) . Since 
x(j) may be of the form cr/7* for integer 1 < 7* < and integer a, we multiply the 
partial solution by 7* to obtain an integral partial solution satisfying 

n— 1 n— 1 

*(n) < 7* n t ^ ri T - 

i=j+l i=j 

Observe that Ay = Ay, so at the end of the procedure we have a solution vector X 
with 



n-1 n-1 / 

< x(l) < ... < x(n) < IJ T< = Xi^df- 1 = ^_ l 



Notice that, since < y(l) < ... < y(n), we also have < x(l) < ... < x(n), 
because of the definition of Ay. ■ 
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In the following we present an example to illustrate the procedure. 



III. An Example 

Let T = {X|AX < 0}, where A is a 4 x 4 matrix and it is only given that its elements 
are in {—1, 0, 1}. In addition, it is known that Y = [2, 3, 7, 29] is an integral solution of 
T. Now, we find an integral vector in T such that its maximum element is 8, which 

satisfies the proposed bound, namely, it is less than 



(2d) 2 



= 16. 

(2=1,71=4 

Recall that T ± — 1,T 2 = 2 and T 3 = 8. We construct a solution [x(l),x(2),x(3),x(A)] 
for A y when y(l) = 2, y(2) = 3, y(3) = 7 and y(4) = 29 . 

1) Initialize x(4) = 1. 

2) A y is the polyhedral cone with constraints c 3 x(3) + c 4 x(4) < 0, where c 3 and c 4 are 
integers with |c 3 | < 8, |c 4 ] < 8, and c 3 y(3) + c 4 y(4) = 7c 3 + 29c 4 < 0. Observe that 
the defining inequalities for A y for this example are < x(3),4x(3) < x(A) and 
5x(3) > x(4); all the other inequalities that we consider are less restrictive. 

3) Initialize x(3) = 1/4. Multiply the solution by 4 to obtain the integral partial 
solution x(3) = 1 and x(4) = 4. 

4) A y is the polyhedral cone with constraints c 2 x(2) +c 3 x(3) +c 4 x(4) < 0, where c 2 , c 3 
and c 4 are integers with |c 2 | < 2, |c 3 | < 2, |c 4 | < 2, and c 2 y(2) + c 3 y(3) + c 4 y(4) = 
3c 2 + 7c 3 + 29c 4 < 0. Observe that the defining inequalities for A y for this example 
are < x(2),2x(2) < x(3) and 2x(2) + 2x(3) < x(4). 

5) Initialize x{2) = 1/2. Multiply the solution by 2 to obtain the integral partial 
solution x{2) = 1, x(3) = 2 and x(4) = 8. 

6) Ay = Ay is the polyhedral cone with constraints cia;(l)+c 2 a;(2) + c 3 a;(3)+c 4 x(4) < 0, 
where ci,c 2 ,c 3 ,c 4 G {-1,0, 1}, and ciy(l) + c 2 y(2) + c 3 y(3) + c 4 y(4) = 2ci + 3c 2 + 
7c 3 + 29c 4 < 0. Observe that the defining inequalities for A y for this example are 
x(l) > 0, x(l) < x(2),x(l) + x(2) < x(3) and x(l) + x(2) + x(3) < x(4). 

7) Choose x(l) = l,x(2) = l,x(3) = 2 and x(4) = 8. 

The output of the algorithm is X = [1, 1, 2, 8]. Thus, X is another integral vector in T. 



IV. Applications 
Some applications of the algorithm are as follows. 

• An important feature of our algorithm is that we do not need to know matrix A to 
produce a bounded solution for the linear program. In many situations, we have 
a solution that satisfies the requirements of a linear system, but we want another 
solution with some upper bounds on its size. It might be impossible (or very time 
consuming) to measure entries of matrix A. Using this method, we can have the 
bounded solution without knowing A. 

• In some streaming applications, we might have the unbounded solution in advance, 
but matrix A arrives later as a query, and we have to produce a bounded solution 
based on A. In streaming problems, we have very limited time for processing the 
query. In our approach, we can solve the problem without the knowledge of query 
in advance. When the query arrives, we already have the answer. 
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• Another interesting problem, that our bounds are useful for it, is the ellipsoid 
algorithm. Assume that we want to find a feasible solution for a linear program. 
We have to start with an ellipsoid that contains at least one feasible solution in the 
core. We then iteratively narrow down, and find a smaller ellipsoid. We have to start 
with an ellipsoid that contains some feasible non-zero integral solution of our linear 
program. In fact there has to be a non-zero integral point in the intersection of our 
polyhedral cone and the starting ellipsoid. If we start with a very small ellipsoid, 
we might not satisfy this property. If the starting ellipsoid is very large, we can 
not prove good bounds on the running time of the ellipsoid algorithm. In this case, 
we have to know some bounds on the size of the starting ellipsoid. Our algorithm 
can be used to get some bounds on the size of the starting ellipsoid. Our current 
bounds do not give polynomial bounds on the running time of ellipsoid algorithm, 
but we hope this approach can lead to such bounds by some polynomials. 

• Another potential application of our algorithm is for finding the routing capacity 
regions of networks. It is known that routing capacity regions of networks can 
be characterized using Farkas lemma as the solution set of infinite set of linear in- 
equalities. But, our algorithm gives an upper bound (finite) on the set of inequalities 
needed to characterize the capacity regions, 
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